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trices, using row and corresponding column permutations. It is a realization of the 
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DESCRIPTION 


Introducti on . This program, REDUCE, reduces the bandwidth and profile 
of sparse symmetric matrices, using row and corresponding column permutations. 
It is a realization of the algorithm described by the authors in [ 4 ) . It 
was extensively tested ind compared with several other programs (5] and was 
found to be considerably taster than the others, superior for bandwidth 
reduction and as satisfactory as any other for profile reduction. 

Outlin e of t he Method . Only an outline of the algorithm is given 
here: a detailed descript 4 can be found in (4J. The algorithm can 

test be described in terms of the adjacen c y graph , C, which has the charac- 
teristic that there is nn edge in C between vertices v^ and v^ if 
and only if a^ i 0 and i 4 j . 

1 . Find the endpoints of a pseudo-d i •’meter; i.e., a pair of 
vertices that are at nearly maximal distance . This is done by a 

finite, iterative process of determining a vertex that is a maximum 
distance away from a given vertex. 

Step 2 . Civen pseudo-diameter endpoints u and v of distance 

k apart, partition the set of vertices into levels L^, L^,...,L^ such 

that adjacent vertices in C are in the same or adjacent levels and 

such that max |L.| is nearly minimized, 
i 1 

S.ep 3 . Number the vertices of G, level by level, beginning at 
an endpoint of the pseudo-diameter. 

Matrix Data Structure . Sparse matrices are typically stored in some 
compact form which takes advantage of the sparsity. The data structure 
assumed here is one which is commonly used in bandwidth and profile schemes; 
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e.g., [11. [2], [31 and [6]. Subroutine REDUCE accepts as input a connect ton 

« 

t able . C, representing the Indices of the nonzero elements of the n x n 
matrix A. The connection table has n rows and m columns where m 
is the number of off-diagonal nonzero elements in the row which has a 
maximum number of off-diagonal nonzero elements (l.e., the maximum degree 
of the graph G). The entries In row 1 of C are the column indices 
of the nonzero elements in row 1 of the matrix A. For example, if 


A - 


X X 0 0 x" 

X X 0 X X 

0 0 X X 0 

0 X X X 0 

X X 0 0 X 


where X represents a nonzero element, then 


C - 


2 5 0 

1 4 5 

4 0 0 

2 3 0 

12 0 
L_ 


The order of indices in a row of C is immaterial. The n-.azero elements 
of the matrix A are never needed, only their indices. 

Test Results . REDUCE was tested on an IBM System/ 360 (model 50) 
computer using the FORTRAN IV G and H compiler? and on a 'DC 6600 computer 
using the FORTRAN (RUN) end FORTRAN extended (FTN) compilers. 

In another paper [51, the authors compared the execution times, band- 
widths and profiles produced by REDUCE with those of five other progra’ns 
on a wide range of problems. REDUCE typically produced the smallest 
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bandwidths; it produced profiles which were on average as small as those 
for any other program; and REDUCE was faster than all the others by at 
least an order of magnitude. 
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ALGORITHM 


c 

c 

c 

C 

c 

* 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

•c 


SUOROUT INC RE DUC C (NOSTK.NR. I OLD . RENUM . NDEG . L VL »L VL S 1 .LVLS2* 

CCSTOR.IBw2.IPF2) 

Subroutine reduce determines a row and coluv*-* permutation which, 

when APPLIED TO a GIVEN SPARSE MA*R|X, PROOUCEo A PERMUTED 
MATRIX WITH A SMALLER BANDWIDTH AND PROFILE. 

The INPUT ARRAY |S A CONNECTION TABLE WHICH REPRESENT? ThE 
INDICES OF THE NONZERO E L E Mr NTS OF THE MATRIX, A. THE ALGO- 
RITHM is DESCRIBED IN TERMS OF THE ADJACENCY GRAPH WHICH 
HAS the CHARACTERISTIC THAT TH*RF IS an EDGE (CONNECTION) 
BETWEEN NODES I AND J IF A(|,J) ,Nf. 0 AND I .NE. J. 

DIMENSIONING INFORVAT » ON--THE FOLLOWING INTEGER ARRAYS MUST BE 
DIMENSIONED IN the CALLING routine. 

NDSTF (NR ,DI ) 

I OLD 02 ) 

RENUM ( 02t I ) 

NDEGID2 ) 

LVL (D2 ) 

L VLSI ( D2 > 

LVLS2 (02 ) 

CCSTOR (02 ) 


D 1 IS .GF. MAXIMUM DEGRFE OF ALL NOOF S . 

02 AND NR ARE .ge. the total NUMBER OF 

NODES IN THF GRAPH, 

STORAGE REQUIREMENTS CAN BE SIGNIFICANTLY 
DECREASED FOR IPM 360 AND 370 COMPUTERS 
by replacing integer nostk by 

I NTEGER* 2 NDSTK IN SUBROUTINES REDuCF, 
OGREE. FNDIAM. tree and NUMBER. 


COMMON INFORMATION — THE FOLLOWING COMMON BLOCK MUST BE IN THE 
CALLING ROUTINE, 

COMMON/GRA/N. IDPTh. I DEG 


EXPLANATION 

NDSTK— 


NR- 

i old ( i ) - 


N- 

IDEG- 


OF INPUT VARIABLES 

CONNECTION TABLE REPRESENTING GRAPH. 

NDSTK ( 1 .J) .NODE NUMBER OF JTh CONNECTION TO NODE 
NUMBER I. A CONNECTION OF A NODE TO ITSELF IS NOT 
LISTED. EXTRA POSITIONS MUST HAVE ZERO FILL. 

ROW DIMENSION ASSIGNEO NDSTK IN CALLING PROGRAM, 
NUMBERING OF I TH NODE UPON IN*=UT, 

IF NO NUMBCRING EXISTS THEN IOLO(I)«|, 

NUMBER OF NODES IN GRAPH (EOuAl TO ORDER OF MATRIX 
MAXIMUM DEGREE of ANY NODE IN THE GRAPH, 




EXPLANATION OF OUTPUT VARIABLES — 

RENUM ( I ) - THE V".W NUMBER FOR THE ITH NCO*.!. 

NDEG(I)- THE DEGREE OF THE ITH NODE. 

I Bw2- THE BAnDwIOTh AFTER RENUMBERING. 

I PF 2- ThE PROF|i_E AFTE* RENUMBERING. 

IDPTH- NUMBER OF LEVELS IN REDUCE LEVEL STRUCTURE. 

THE FOLLOWING ONLY HAVE MEANING if The graph WAS CONNECTED-- 
L VL ( I ) - INDEX INTO LVLS1 TO THE F I RSI NODE IN LEVEL 1. 

LVL ( I ♦ I >-L VL ( I ) * NUMBER OF NODES IN |TM LEVEL 
L VL SI— NODE NUMBERS LISTED BY LEVEL. 
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C LVLS?(I>“ The level assigned to nooe I pv reducc* . 

c 

C FORKING STOWAGE VARJABlE — 

C CCSTOR 
f 

C local STORAGE — 

C COMMON/CC/— SUBROUT I NE S REDUCE. SORT2 AND PIKLVL ASSUME THAT 

c the graph has at most so connected components, 

C SUBROUTINE E NO I AM ASSUMES THAT THERE ARE AT MOST 

C 100 NODES IN THE L a <T i.FVFL, 

C COMMON/LVLM/-SUI3ROUT I NCS SETUP AND PIKLVu ASSUME THAT THERE 
C ARE AT MOST 100 LEVELS, 

C 

c 

C USE INTEGER«2 NDSTK w|TH An I BM 360 OR 370, 

c 

INTEGER NrSTK 

INTEGER STNODE .RVNODE . WENUM . xC . SORT 2 . STNUM . CCST OR • S I 2E . STPT .SP.NUM 
COMMON /gra/ n,igrth*ideg 

IT IS ASSUMED THAT THE GRAPH HAS AT MOST 50 CONNECTED COMPONENTS, 

COMMON /CC/ XC. SIZE <50 » .STPT (50) 

DIMENSION CCST OR ( I ) . I OLD ( 1 ) 

01 MENS ION NDSTK (NR. I i • L VL < I ) .LVLS1 (I ) » LVLS2 ( I ) .RENUM < | ) .NDEG < I I 
I Bw2*0 
iPr2»o 

SET RENi/M ( I > >0 FOR ALL I TO INDICATE NODE I IS UNNUMBERED 
DO I O I * I .N 
RPNUM ( f ) «0 
10 CONTINUE 

COMPUTE DEGREE O- EACH NODE AND ORIGINAL BANDWIDTH AND PROFILE 
C 

call DGREE (NDSTK.NR.MDEG. IOLD. IBWI . I PF | ) 

C SBN'JMb LOW END OF AVAILABLE NUMBERS FOR RENUMBERING 
C STNUM* HIGH END OF AVAIlAE 3;F NUMBERS FOR RFNUMdEW I NG 
SBNUM* l 
S t NUM=N 

C NUMBER THE NODES OF DEGREE ZERO 
DO 40 I *1 «N 

IF (NDEG < I ) «GT .0 ) GO TO «0 
RENUM ( | ) «STNUM 
STNUM«STN* JM- i 
40 CONTINUE 

C FIND AN UNNUMBERED NODE OF M|N DEGREE TO START ON 
50 LOwDG* I DEG+ 1 
NFLG« 1 
I SO I R* I 
DO 70 I » I .N 

IF (NDEG ( I ) ,GE «lOwOG I GO TO 70 
IF (RENUM ( | >,GT,0) GO TO 70 
LOWDG-NDEG ( I ) 

STNODE* I 
70 CONTINUE 
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t Ffvo psfudo— o i a me rER and associated level strictures. 

C ST?, -DC AND R VNODE ARE THE ENOS OF THE DIAM AND LVLSl ANO LVLS2 
C ARE THE RESPECTIVE LEVEL STRUCTURES. 

CALL END I AM ( S TNOOE • RVNOOt • HOSTS • NR « NDCG »L VL »L VLS 1 .1 VLS2.CCSTOR. 

- IDELT ) 

I E < NOE G I ST NODE ) *LE • NOE C ( RVNOOE I I GO TO 75 
c nelg indicates the end to begin numbering on 

NELG»-1 
STNOOL =RVNOOE 

75 CALL SETuP<LVL .L VLSI .L VLS? ) 

C FIND all The CONNECTED COMPONENTS <XC COUNTS THEM) 

XC*0 

lroot«i 

LVLN» I 

DO 00 | • I, N • 

IE(LVL< I t.NE.O) GO TO BO 
XC«XC* I 

STPT ( XC ) ^LROOT 

CALL TREE I I .NOETK*NR.LVL»CCSIOR.NDLG.LV I _wT» .LV'.BOT .LVLN -.MAXLW .N ) 
SIZE < XL l «LVLBOT4LVL»TH-LROOT 
LROOT «LVLBOT4LVLWTH 

lvln*lroot 
BO CONTINUE 

IE CS0RT2(DMY l.CO.O) GO TC 90 

CALL PIELVLlLVLSl .LVLS2.CCST0R. I DEL T . I SO I R ) 

C ON RETURN FROM PJXLVL. I SD I R INDICATES THE DIRECTION THE LARGEST 
C COMPONENT FELL. I SD I R IS MODIFIED NOW TO I NO 'CATE THE NUMBERING 
C DIRECTION. N'JM IS SET TO 1 HE PROPER VALUE FOR THIS DIRECTION. 

90 I SDIR= I SD|R*NELG 
NUM.SBNUM 

IF ( ISDIR.LT.O ) NUMs STNUM 

CALL NUMBER ( STNODE.NuM.NDSTK .LVLS2.NDEG.RENUM.LVLSI .LVL .nr.nelg. 
-IB*?. IPE2.CCSTOR. I SO I R ) 

C UPDATE STNUM op SBNUM after NUMBERING 
1F( ISDIR.LT.O) STNlHsNuM 
|E( ISCIR.GT.O) SBNuMsNuM 
IE (S0NUM.LE. STNUM ) GO TO 50 
IF C IBV2.LE. IBWI ) RETURN 

C IE ORIGINAL NUMBERING IS BETTER THAN NEW ONE. SET v»P TO RETURN IT 
DO 100 I « I .N 

RENL'M < | )• | OLD < I ) 

100 CONTINUE 
!BW2> IBW! 

|PF?*|PFI 

RETURN 

END 
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SUBROUTINE dgqee cndstk.nr.noeg. iolp. ibwi • ;pfi 

c 

C DOPE E COMPOTES The DEGREE OF EACH NODE IN NDST* AND STORES 
C IT IN The ARRAY NDfcG. T Hfc BANDWIDTH AND profile fop The original 
C OP INPUT RENUMBERING OF THE G»APM IS COMPUTED ALSO. 

c 

C USE INTEGER*? NDST< WITH AN IBM 360 op 370. 

C 

INTEGER NOSTr. 

COMMON /GRA/ N IDPTh.IDEG 

DIMENSION ND3T< CNR, I ) . NOEGS I > . IC.DC I > 

I Bwl *0 
|PF| nO 

DO ICO l»i.N 
NOEGC I ) «0 
I Rw ■ 0 

DO 80 .1*1 . I DEG 
| tst*nostf t I . J ) 

IF C I TST ) 90.90.50 
50 NDFGC |)*NOEGC I )4 1 

ID IF * I OLD C I I - I OLD < I T5T | 

IF ( IRW.LT. IDIF | |Rw* I D I K 
80 CONTINUE 
90 |PE | a |PF| 4|Rw 

irclRw.G~.IBwl) I Bw I * I &W 
IOO CONTINUE 
rf turn 
END 
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SUBROUT INE FNOIAMISNDl • SND2 . NDSTK .NR .NOCG .LVL • LVL S 1 »L VLS2 » 

-I WK. 10FLT ) 

c 

C FNDIAM is THl CONTROL PROCEDURE FOR FINDING THE PSEUDO— D I AWE TER OF 
C NDSTK AS WELL AS THE LEVEL STRUCTURE FROM E*CH END 
C 

C SND I - 

c 
c 

C SN02- 
C LVLSl - 
C LVLS2- 
C IPFLT- 

c 

C LVL . I W< 

C 

C USE INT~SER*2 NDSTK *|Th AN IBM 3€>0 OR 370. 

C 

INTFGFP NDSTK 

INTEGER FLAG.SND.SND1.SND2 
common /gra N.IDPTH 

IT IS ASSUMED THAT THE LAST LEVEL HAS AT MOST 100 NODES. 

COMMON /CC/ NOLST(IOO) 

D I MENS I ON NDSTK <NR. 1 | . NOFG < I ) »L L I | ) .LVLSl ( I J .LVLS2 ( I ) • | WK < | ) 

flag*o 

M T W2*N 
SND* SND 1 

C ZERO LVL TO INDICATE ALL NODES ARE AVAILAOLF TO TREE 
20 DO 25 ! * l « N 
LVL ( I »«0 
25 CONTINUE 
L VLN* 1 

C DROP A TREE FROM SND 

CALL TREE (SND.NOSTK .NR. LVL » I WK . NDEG »L VL WT H,L’'LbOT . L VLN. M AXL W . MT W2 ) 
IF (FLAG.GE.l ) GO TO 110 

flag*i 

70 I DPI H*L VLN— 1 
MTwl =MAXLW 

C COPY LFVEL STRUCTURE INTO LVLSl 
DO 75 I * l *N 

LVLSl < I ) - LVL ( I ) 

75 CONTINUE 
NDXN=1 
NOXL *0 
MTW2*N 


CN INPUT THIS IS THE NODE NUMBER OF THE FIRST 
ATTEMPT AT FINDING A DIAMETER. ON OUTPUT IT 
CONTAINS THE ACTUAL NUMBER USED. 

ON OUTPUT CON i A I NS OTHER END OF DIAMETER 

ARRA ( CONTAINING LEVEL STRUCTURE WITH SNOl AS ROOT 

ARRAY CONTAINING LEVEL STRUCTURE WITH SND2 AS ROOT 

flag used in picking finAc. level structure, set 

*1 IF WIDTH of LVLSl .LE. WIDTH OF LVLS2 . OTHERWISE *2 
WORKING STORAGE 





C SOPT LAST LEVEL by DEGREE AMD STOWE IN NDL£T 

call SO«TDG<Nl>LSr. I WK (LVLBOT ) .NGXL.LVLWTM.NDEG » 
SNO»NDLST < I ) 

GO TO 20 

110 IE ( IDPTM.GE.LVLN- I > GO TO 120 
C START AGAIN WITH NEW STARTING NODE 
SNOI «SND 
GO TO 70 

• 20 IE (MAXLW.GE.HTW2 ) GO TO 130 
MTW2» VAxL W 
SND2*SN ' 

C STORE NARROWEST REVERSE LEVEL STRUCTURE IN t VLS2 
T>0 I ?•> I = 1 • N 

LVLS2 ( I I «LVL ( I ) 

12* '•ONTINUE 

•Sw IE (NDXN.EQ.NDXL # GO TO 140 

c try next node in ndlst 

NOXN*NDXN+ 1 
SNC = NDL ST ( NP XN ) 

GO TG 20 
1 AO I DrLT* I 

IE (MTW2.LE.MTW1 ) IDELT*2 

RETURN 

ENO 
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SUBROOT INE TREE I I ROOT • NDSTK • NR .L VL » I WK .NDEG.LVLWTH.LVLBOT » 
-LVLN.MAXLW. IBOPT > 

c 

C tree DROPS » TREE IN NDSTK FROM I ROOT 

c 

C L VL “ 

C 

c 

C I WK- 
C 
c 

C LVLWTH- 
C LVLBOT- 

c 

C MAXLW- 
C LVLN- 
C 

c 
c 

C 1BORT- 


iNTEGER*2 NDSTK WITH AN IBM 360 OR 370. 
integer NOFTK 

DIMENSION NDSTK { NR . 1 ).LVL< I ) • i WK ( 1 ) *NDEG< 1 ) 
M A XL W* 0 
I TOPtLVLN 
I NOW S L VLN 
lvlbot^lvln 

LVLTOP*LVLN4- I 
L VLN* I 

L VL ( I ROOT >*1 
IWKf I TOP)* I ROOT 
L VL N=LVLN* 1 
I .NOW: 1 w< ( I NOW ) 

NDROW»NDEG ( I WK NOW ) 

DO 40 J*1 .NOROW 

! TEST =NDSTK ( | WKNOW • J ) 

IF (LVL I I TEST ) .NE.O ) GO TO 40 
LVL I I TEST J *LVLN 
I TOP= I TOP* I 
I WK ( I TOP ) * I TEST 
CONT INUE 
I NOW* I NOW* I 

I F ( I NOW .LT • LVL TOP ) GO TO 35 

LVLWTH=LVLTOP-LVLBOT 

IF <I4AXLW*LT*UVLWTH) maxlw*lvlwth 

IF (MAXLW.GE. IBORT ) RETURN 

IF C I TOP.LT. LVL TOP ) RETURN 

LVLBOTs I NOW 

LVLTOP: ITOP^ I 

GO TO 30 

ENO 


C 

r 

C USE 
C 


30 

35 


40 


ARRAY INDICATING AVAILABLE NODES IN NDSTK WITH ZERO 

ENTRIES, tree FNTFRS LEVEL NUMBERS ASSIGNED 

DURING FXFCUTICN OF Th|S PROCF0URF 

ON OUTPUT CONTAINS NODE NUMBERS USED IN TREE 

ARRANGED BY LEVELS * I WK (L VLN ) CONTAINS I ROOT 

AND I WK (L Vl.BOT 4-L VL WTH- 1 ) CONTAINS LAST NODE ENTERED) 

ON OUTPUT CONTAINS WIDTH OF LAST LEVEL 
ON OUTPUT CONTAINS INDEX INTO I WK OF FIRST 
NODE IN LAST LEVCL 

ON OUTPUT CONTAINS THE MAXIMUM LEVEL WIDTH 

ON INPU* THE FIRST AVAILABLE LOCATION IN I WK 

USUALLY ONE BUT IF IWK is USED TO STORE PREVIOUS 

CONNECTED COMPONENTS. LVLN IS NEXT AVAIL AbLE LOCATION. 

ON CUTPUI THE TOTAL NUMBER OF LFVELS ♦ I 

INPUT PARAM WHICH TRIGGERS early RETURN IF 

MAXLW BECOMES .GE. IBORT 
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SUBROUT INF SORTOGISTKl . STK2.X1 • X2 «NDEG ) 

c 

C SORTOG SORTS STK2 BY DEGREE OF THE NODE AND AOOS 
C OF STK1 IN ORDER OF LOWEST TO HIGHEST DEGREE* XI 
C NUMBER OF NODES IN STK I AND STK2 RESPECTIVELY. 

C 

INTEGER XI «X2.STKI .STK2. TEMP 
COMMON /GRA/ N.IDPTH 
D I MENS I ON NDFG < I ) .STXl ( I > .STK2< 1 ) 

I ND r X2 
10 I TEST* 0 
I NO* IND-1 

IF I INO.LT. 1 » GO TO *0 
DO 30 I * l « I NO 
J« 141 

|STX2»ST<2 ( I 1 
JS T <?*STK2 ( J ) 

IF<NDEG< ISTK 2 ).LE*NDEG(JSTK2) ) GO TO 30 
I TEST * I 
TFMP«ST<2( I 1 
STK2I I ) *STK2 ( J ) 

STK2 ( J) sTEML 
30 CONTINUE 

IF ( ITEST.EO* 1 1 GO TO 10 
4 0 DO SO I ■ I « X2 
XI *K1 +1 

STXl (XI )«STX?( I ) 

50 '•ONTINUE 
RETURN 
ENO 


IT TO THE ENO 

and x2 are the 
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SU6RCUT INE SETU' 3 <LVL .LVlSI .L VLS2 ) 


c 

C SETUP COMPUTES THE REVERSE LFVF'.ING INFO FROM LVLS2 A NO STORES 
C IT INTO LVLS2. NACUM(I) IS INITIALIZED TO NODFS/ITm LEVEL FOR NOOES 
C ON TMF PSEUOO-DJAMETER of THE GRAPH. LVU IS INITIALIZED to NON- 
C ZERO FOR NOCES ON THE PSEUOO-OIAM ANO NO! IN A DIFFERENT 
C COMPONENT OF THE GRAPH. 

C 

COMMON /GRA/ N.IDPTH 

IT IS ASSUMED THAT THERE ARE AT MOST 100 LEVELS. 

COMMON /LVLW/ NH| GH( 1 00 ) <NLOW( I OG ) t NACUMt | 00 ) 

D I MENS ION L VL ( I I .L VL S| (1 ) . L VL S2 f I ) 

DO 30 I *1 . IOPTh 
NACUM ( | ) = 0 

30 CONTINUE 

DO 1 *0 I ■ I . N 
L VL I I ) * 1 

LVLS2 ( I I * IDPTH4 | -LVLS2 ( I ) 

1TEMP=LVLS2( I ) 

IFUTEMP.GT. IDPTH) GO TO 140 
IF I ITEMP.NE.LVLS1 ( I > I GO TO 100 
NACUM ( | TEMP ) sNACUM | | TEMP )♦ I 
GO TO l«0 
I 00 LVL ( I > *0 
140 CONTINUE 
RE TURN 
END 
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INTEGER FUNCTION S0PT2I0MV) 


c 

C S0PT2 SORTS SIZE A NO STPT INTO DESCENDING ORDER ACCORDING TO 
C VACUCS OF SIZE. XC -NUMBER OF t'NTPIES IN EACH ARRAY 

c 

INTEGEP TEMP. ccstop. size. stpt.xc 

IT IS ASSUMED THAT the C.RAPm HAS AT most SO CONNECTED COMPONFNTS. 

COMMON /cc/ XC .SIZE* SO J » STPT I SO ) 

SOPT2=0 

IF (XC.EG.O ) PE TURN 
SORT Z 8 1 
I NO * X C 
10 I TEST * 0 
I NO« I NO- I 

I r I ; \ . L T . I > PC TURN 
DO I 7 , S | . I NO 
j«m 

IF (SIZE C I I.GE.SIZF (Jl ) GO TO 17 
I TEST * I 
TpMP= SIZE ( I ) 

S I ZE ( I I *S I ZE ( J ) 

SIZE ( J) »TEMP 
TEMP. STPT < I ) 

STPT ( I ) -STPT ( J | 

STPT ( J j «TEMP 
17 CONTINUE 

IE ( ITEST.EO. 1 ) GO TO 10 
PE TURN 
ENO 
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subrout ine pixlvl clv-.si .lvlS2*ccstor* idflt, isdib > 

t 

c 

C PlXLVL CHOOSES the LEVEL STRUCTURE USED IN NUMBERING GRAPH 
C 

C LVLSI- 
C LVLS2- 

c 

C CC5.TOR- 
C ICfLT- 
C NH | GH 
C NLOW- 
C NACUM- 
C XC- 
C SIZE <11- 
C ST°T C I > — 

C 1SDIR- 

C 

INTEGER CCSTOR.SIZE«STPT,xC.END 
COMMON /GRA/ N.IDPTh 

c 

C IT IS ASSUMED THAT THE GRAPH HAS AT MOST *."> COMPONENTS AND 
C THAT There ape at most 100 levels* 
c 

COMMON ✓LVl-W/’ NH | GK ( too ) .NLOW < I 00 1 .NACUM ( 1 00 ) 

COMMON /cc/ XC «S1 ZE < 50 ) .STPT (50 ) 

DIMENSION LVLSl M ) .L VLS2 C 1 ) .CCSTORl I 1 
C EOR EACH CONNECTED COMPONENT DO 
DO 270 1*1 • XC 

j.T.TPT ( I ) 

END* S I ZE < I 14J-1 

C SET MHIGH AND NLOW EQUAL TO NACUM 
DO 205 X* l , 1DPTH 
NHIGMIX) *NAClJM(X ) 

NLOW < X ) ^ NACUM (X ) 

205 CONTINUE 

C update NHIGH AND NLOW FOR EACH NODE IN connected COMPONENT 
DO 21 0 XsJiEND 
1 NODE = CC STOB ( X ) 

LVLNH=LVLSl < INODE ) 

NH | GH (LVLNH ) * NH | GHILVLNH )♦ 1 
LVLNL *LVLS2 C I NOOE ) 

NLOw ILVLl*:L > *NLO* CLVLNL ?♦! 

210 CONTINUE 
MA X 1 * 0 
MAX2*0 

C SET MAX1*LARGEST NEW NUMBER IN NH I GH 
C SET MAX2*LABGEST NEW NUMBER IN H.OW 
DO 240 X« 1 • IDPTH 

IE < 2*NACUM (X ) .EO.NLOW (X )^NH| GH(X ) ) GO TO 240 
!F(NHiGH(K).GT,MAXi 1 MAX 1 =NM 1 GH C X > 

IE (NLOWIX » *GT*MAx2 ) MAx2*NLOW(K> 

24C CONTINUE 


on input contains forward leveling info 
on input contains reverse leveling ineo 
ON OUTPUT the final LEVEL STRUCTURE chosen 
on input contains connected COMPONENT info 

ON INPUT *1 IE WDTH LVLSI»wDTH LVLS2. «2 OTHERWISE 

XEEPS TRACX or LEVtL WIDTHS FOR H I GH NUMBERING 

XtEPS TPACX OF LEVEL WIDTHS FOR LCw NUM«tS|NG 

XEEPS 1RACX OF LEVEL WIDTHS FOR CHOSrN LFVEL STRUCTURE 

NUMBER OF CCNNFCTFD COMPONENTS 

SIZE OF |TH CONNCCTEO COMPONENT 

INDEX INTO CCSTORE OF 1ST NODE IN I TH CON COMPT 
flag which indicates WHICH WAV ThE LARGEST connectfo 
COMPONENT FELL* * ♦ I IF LOW AND -I IF HIGH 


J.5 


C SET IT« NUMBER O^ LEVEl. STRUCTURE TO BE USED 
IT«I 

|F»MAXI .GT.MAX2 > I T B 2 
IF (MAXI «EQ»MAX2 > IT«IDFlT 
I F I IT.E0.2 ) GO TO 265 
IF < l .EG. 1 » I SO I R*— I 

c copy LVLSI into lvlS2 for each node, in connected component 
DO 260 (C* J «END 
INODE =CCSTOR<< ) 

L VLS2 ( I NODE > «L VL S I < i NODE ) 

260 CONTINUE 

c update nacum to be the same as nhigh 

OO 262 K * I • IDPTH 

NACU»» «K I : NM I fjM ( K ) 

262 CONTIN'JE 
GO to 270 

C MDATE NACUM TO BE THE SAME AS NLOV 
265 DO 267 E*1. IDPTH 

NAC.UM <K ) »N| OW (K I 
267 CONTINUE 
2 70 CONTINUE 
RETURN 
END 
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SUBROUT IKE NUMBER ( SND • NUM , NDSTK .LVLS2»ND€G.RENUM.LVlST .LSTPT .MR# 
-NFLG. IBW2. IPF2. 1PFA. ISOIR) 

NUMBER PRODUCES THE NUMBER INC, OF THE GRAPH FOR MIN BANDWIDTH 
C 

C SNO- ON INPUT THE NOOF TO BEGIN NUMBERING ON 

C NUM- ON INPUT AND OUTPUT. T Hr NEKT AvAlLABLF NU m HER 

C LVLS2- THE LEVEL STRUCTURE TO BE USED IN NUMBERING 

C RFNUM- T He ARRAY uSF D TO STORE THE ME * NUMBERING 

C LVLST- ON OUTPUT CONTAINS LEVEL STRUCTURE 

C LSTPT ( | * - ON OUTPUT. » NDEX INTO LVlST TO FIRST NODE IN ITH LVL 

C LSTPT (!♦!» •• LSTPT' I ) * NUMBER OF NODES IN ITH LVL 

C NFLG- = ♦ I IF SND S FORWARD LNO OF PSFuDO-DIAM 

C *-l IF SND IS RFVERSE END OF PSFuDO-DIAM 

C IBw2- BANDWIDTH OF NE W NUMBERING COMPUTED By NUMBER 

C IPE2- PROFILE OF NEW NUMBERING COMPUTED BY NUMBER 

C I PF A— . ORK | NG STO-JAGi* UjED TO COMPuTF PROF |t_E AND BANDWIDTH 

C ISDIR- INDICATES STEP DIRECTION USED IN NUMBFRINGC+1 OR -l | 

C 

C USF INTEGEO»2 NDSTK WITH an IOM 360 OR 3?0* 

c 

INTEGER NDSTK 

INTEGER SND. STFA.STKB. ST«C. STKD. XA.XB»XC.XO«CX .ENO.RENUM .TEST 
COMMON /GBA/ N.IDPTH.IDEG 
C 

C the STORAGE IN COMMON OLOCKS CC AND LVL* IS NOW FREE AND CAN 
C BE USED FOR STACKS. 

C 

COMMON /LVL W/ STKA ( 1 OO ) , STXB ( l 00 ) .STXC ( I 00 ) 

COMMON /CC/ STKD(IOO) 

DIMENSION IPFA ( 1 ) 

D I MENS I ON NDSTX ' NR . I ) .LVt.S? ( I » » NDEG I I ) .RENuM ( 1 ) .LVLST ( I ) .LSTPT! I ) 
C SET UP LVLST AND LSTPT FROM LVLS2 
DO 3 I » I • N 
I PEA (||*0 
3 COwT I NUE 
NSTPT* I 

DC B 1*1. IOPTH 
LSTPT I 1 ) »NSTPT 
DO S J*1 N 

IFILVLS. (J).NE.I > GO TO 5 
LVLST (NSTPT )* J 
NSTPTsNSTPT-M 
5 CONTINUE 

LSTPT ( I DPTH+ I »*NSTPT 
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C STKA. STKB. STKC A NO STkD ARF STACKS WITH POINTERS 
C XA.XR.XC. AND XD. CX IS A SPFCIA*. POINTER INTO STKC WKJCH 
C INDICATES THE PART I CUk_ AW NODE BE I NG PROCESSED. 

C LVlN KEEPS TRACK OF TmE LEVEL wE ARE WORKING AT. 

c initially stkc contains only the initial node. snd. 

LVLN»0 

IF INFLG.l T.O ) lvln*ioptm*i 
XC»I 

STKC I XC l«SND 
10 CX»| 
xr>» o 

lvln«l vlnknflg 

LST«LSTPT (LVLN ) 

L ND*L J ip r ILVLN* I >-l 
C BEGIN PROCESSING NOOE STKC (CXI 
20 |PRO«STKC (CX I 
RENU*M I PRO I «NUM 
NUM = NUM* I SO ! R 
E»..>«NOFG( POO) 

XA«0 

XB«0 

C CHECK ALL ADJACENT NODES 
DO SO | » I ,ENO 

TcT.T»NDSTK ( I PRO. I I 
| NXsPFNUM ( TEST | 

C ONLY NODES NOT NUMRF»ED OR ALREADY ON A STACK ARF ADDED 
IF(INX.EO.O) GO TO 30 
IF(INX.LT.O) GO TO 50 

C DO PRELIMINARY BANDWIDTH and profile calculations 
NR W* (RENUM ( I PRO I- I NX » * I SD I R 
IF ( I SDIR.GT .0 I INxiBENiM ( I PRO I 
IF ( IPFA ( INX ) .LT.NBW I IPFAt |NX| =NHW 
GO TO 50 

30 RFNUM (TEST )»— I 

C PUT NODES ON SAME LEVEL ON STKA « ALL OTHERS ON STK3 
IF ( LVLS2 (TEST I *EO«L VLS2 ( I PRO I I GO TO 40 
XB«XB-M 

STkB< XB > «TEST 
GO TO 50 
40 XA = X A •* l 

STKA ( XA ) * TEST 
50 CONTINUE 

C SORT STKA AND STKB INTO INCREASING DEGREE AND ADO STKA TO STKC 
C AND STKB TO stkd 

IF (XA.EO.O ) GO TO 55 
IF (XA.EO* I ) GO TO 52 
CALL SORTDG(STKC.STKA.XC.XA.NOEG) 

GO TO 55 
52 XC«XC41 

STKC (XC I « STKA (XA ) 

55 I F ( xB.EO. 0 I GO TO 6* 

I F ( XB.EO. I > GO TO 62 

CALL S0RTDG(STKD.STKB.XD.XB.N0EG» 

GO TO 65 
62 XP»XD*i 

STKD(XD)*STkB(XB) 
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C BF SURF TO PROCESS ALL NOOFK JN STKC , 

65 CK*rx+l 

If<xc.gc*cx> no to ?o 

C WHEN STkC IS EXHAUSTED LOOK FO« MlN DEGREE NODE IN SAME LEVEL 
C WHICH HAS NOT BEEN PROCESSED 
MAX* IDI &♦ I 
SND«N^ I 

OO 70 1-LST.LND 
TEST«LVLST ( I | 

|F (RENOM( TEST >*N£.0 ) GO TO 70 
IF (NOEG ( TEST ) «GE *MAX > GO TO 70 
R^NU* ( SN"> ) *0 

PFNUM < TFST I •- I 
KA*»N')Ff.lTf ST > 

SND* TEST 
70 CONTINUE 

|F (SND»EO«N* I ) GO TO 75 
XC- C ♦ I 
STKC ( XC )»SND 
GO T'J 20 

C IF STKO IS fcMPiy wE ARE DONE. OTHERWISE COPT STKD ONTO STKC 
C AND BEGIN PROCESSING NFw STKC 
75 |F(XD*FO. GO TO 1 00 
OO BO I «l .XO 

S T K C I I » »STKDI I I 
60 CONTINUE 
XC«XD 
GO TO 10 

C DO FINAL BANDWIDTH ANO PROrjLE CALCULATIONS 
100 DO 1?0 I ■ I «N 

|F l IPFA ( | ) ,GT. IBw2 > I bw? * I PF A < | ) 

|PF2« I PF 2 ♦ I Pf A ( | ) 

120 CONTINUE 
RFTuRN 
END 
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